Anisotropic Homogeneous Turbulence: hierarchy and intermittency of scaling 

exponents in the anisotropic sectors. 
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We present the first measurements of anisotropic statistical fluctuations in perfectly homogeneous 
turbulent flows. We address both problems of intermittency in anisotropic sectors and hierarchical 
ordering of anisotropies on a direct numerical simulation of a three dimensional random Kolmogorov 
flow. We achieved an homogeneous and anisotropic statistical ensemble by randomly shifting the 
forcing phases. We observe high intermittency as a function of the order of the velocity correlation 
within each fixed anisotropic sector and a hierarchical organization of scaling exponents at fixed 
order of the velocity correlation at changing the anisotropic sector. 
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At the basis of Kolmogorov 1941 theory there is the idea 
of restoring of universality and isotropy at small scales in 
turbulent flows. Memory of large scale anisotropic forc- 
ing and/or boundary conditions should be quickly lost 
during the process of energy transfer toward small scales. 
The overall result being a local recovering of isotropy and 
universality for turbulent fluctuations at small enough 
scales and large enough Reynolds numbers. 

In recent years, a quantitative investigation of restor- 
ing of isotropy in experimental anisotropic turbulence 
jl],D, numerical homogeneous shear flows (^^] and nu- 
merical channel flows || questioned the main Kol- 
mogorov paradigm, speaking explicitly of persistence of 
anisotropies. Some theoretical work has also been done 
Q in order to understand how to properly link the in- 
variance under rotation (SO (3) symmetry group) of the 
Navier-Stokes equations and the analysis of anisotropic 
fluctuations of velocity turbulence correlations. 
The observed anisotropic effects in small scales turbu- 
lence is both a theoretical challenge and a very actual 
practical problem, opening the question whether any, re- 
alistic, anisotropic turbulent flows can ever possess statis- 
tical features independent of the (anisotropic) boundary 
and forcing effects. This goes under the name of univer- 
sality. 

Neglected anisotropic effects in high Reynolds numbers 
flows have also been proposed to be at the origin of 
different statistical properties measured for transversal 
and longitudinal velocity fluctuations [Q. Importance of 
properly disentangling isotropic and anisotropic fluctua- 
tions has also been demonstrated in the analysis of inter- 
mittency in channel flow turbulence j8|. 
Important step forward in the analysis of anisotropic 
fluctuations have recently been done in Kraichnan mod- 
els, i.e. passive scalars/vectors advected by isotropic, 
Gaussian and white- in-time velocity field with large scale 
anisotropic forcing |p[-pd[. 

In those models, anomalous scaling arises as the results 
of a non-trivial null-space structure of the advecting op- 



erator. In these cases, correlation functions in differ- 
ent sectors of the rotational group show different scal- 
ing properties. Scaling exponents are universal: they do 
not depend on the actual value of forcing and boundary 
conditions, and they are fully characterized by the or- 
der of the anisotropy. Non universal effects are felt only 
in coefficients multiplying the power laws. Coefficients 
are fixed, in principle, by requiring matching with non- 
universal boundaries conditions in the large-scale region. 
Similar problems, like the very existence of scaling laws 
in the anisotropic sectors and, if any, what are the val- 
ues of the scaling exponents and what is the dependency 
from universal/non- universal effects are at the forefront 
of experimental, numerical and theoretical research in 
true turbulent flows. Only a few indirect experimental 
investigation of scaling in different sectors Jil|[l6| and di- 
rect decomposition in Channel flow simulations p|,p7|,p[ 
have, at the moment, been attempted. 
The situations is still unclear, evidences of a clear im- 
proving of scaling laws by isolating the isotropic sector 
have been reported, supporting the idea that the unde- 
composed correlations are strongly affected by the su- 
perposition of isotropic and anisotropic fluctuations 
Preliminary evidences of the existence of a scaling law 
also in the sectors with total angular momentum j = 2 
have been reported |l5|,[l6|]. The value of the exponent 
for the second order correlation function being close to 
the dimensional estimates £^ -2 = 4/3, |l8|| . 
All these preliminaries investigation in real turbulent 
flows are flawed by the contemporary presence of 
anisotropies and strong non- homogeneities. The very ex- 
istence of scaling laws in presence of strong non homo- 
geneous effects can be doubted. SO(3) decomposition 
becomes soon intractable as soon as non-homogeneous 
effects cannot be neglected Q. Moreover, in many exper- 
imental situation, anisotropies are introduced by a shear 
forcing coupled to all turbulent scales: something which 
prevents the possibility to study "pure" inertial physics. 
To overcome this problems we performed the first nu- 
merical investigation of a turbulent flow with strong 
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anisotropic forcing confined to large scales and perfectly 
homogeneous on a numerical resolution 128 3 and 256 3 . 
We studied a fully periodic Kolmogorov flow with ran- 
dom, delta correlated in time, forcing phases, which we 
decide to call a "Random-Kolmogorov Flow" (RKF). 
In this letter we present direct measurement of scaling ex- 
ponents in sectors up to total angular momentum j = 6. 
Our main results support the existence of a hierarchi- 
cal organization of exponents, i.e. continuous increase 
of exponents as a function of j. We also found a much 
stronger intermittency in the anisotropic sectors than in 
the isotropic one. We conclude with a few comments and 
proposals for further work in the field. Let us begin to 
expose a few technical details on the simulations. We 
performed a direct numerical simulation of a fully peri- 
odic flow with anisotropic large scales forcing. In details, 
we have chosen a random forcing pointing only in one 
direction, the z axis, with spatial dependency on the x 
direction only on two wavenumbers fci = (1,0, 0),fc2 — 
(2,0,0). Namely: fi(k 1:2 ) = ^,3/1,2 exp^fli^) where 
fi > /2 are two constant amplitude and Q\ , 62 are two 
random phases, delta correlated in time. The random 
phases allows for a homogeneous statistics also in the 
-otherwise- non-homogeneous direction spanned by the 
two wavenumbers, i.e. we have instantaneously a large 
scale non-homogeneity in the x direction which is aver- 
aged out by the time evolution thanks to the random 
re-shuffling of the forcing phases. We studied the RKF 
at resolution 128 3 and 256 3 , we collected up to 200 eddy 
turn over times for the smallest resolution and up to 50 
eddy turn over times for the largest resolution. Such a 
long averaging is necessary because as in any strongly 
anisotropic flow we observe the formation of persistent 
large scale structures inducing strong oscillation of the 
mean energy in time. 

In Fig. 1 we show, for example, a typical time evolu- 
tion for the total energy and total energy dissipation in 
our runs, it is interesting to notice how the high frequen- 
cies oscillation at large scales (total energy) induces by 
the random forcing are completely absent at small scales 
(energy dissipation). 

In order to increase the scaling range extension wc have 
used an hyperviscosity with a squared laplacian. Inset of 
Fig. 1 quantifies our degree of homogeneity. We have a 
high degree of homogeneity (more than 95%) in the two 
transverse directions, y, z, while we still observe small os- 
cillations in the x directions (of the order of 10%); these 
oscillations are due to statistical fluctuations, they must 
averaged out in the limit of infinite statistics. 
Let us now discuss the SO (3) decomposition of longitu- 
dinal structure functions: 

S p (r) = ([(v(x)-v(x + r))-r] p ), (1) 

where we have kept only the dependency on r neglect- 
ing the small non homogeneous fluctuations. We expect 



that the undecomposed structure functions are not the 
"scaling" bricks in the theory. Theoretical and numer- 
ical analysis showed that one must first decompose the 
structure functions on the irreducible representations of 
the rotational group and than asking about the scaling 
behavior of the projection. In practice, being the longi- 
tudinal structure functions scalar objects, their decom- 
positions reduces to the projections on the spherical har- 
monics: 

00 j 
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Where we have used the indeces j, m to label the to- 
tal angular momentum and its projection on a reference 
axis, say z, respectively. The whole physical informa- 
tion is hidden in the coefficients >S^ m (|r|). In partic- 
ular, the main question we want to address here con- 
cerns their scaling properties: S J p m {\r\) ~ Aj m \r\^^ 
and (in the case) what one can say about the values of 
the scaling exponents, and on their robustness against 
large scale physics (universality issue). Theoretical ar- 
guments suggest that if scaling exponents exist they de- 
pend only on the j eigenvalue [|9|. If true turbulence 
follows the Kraichnan models behavior, we should ex- 
pect universality of the scaling exponents (independence 
of large-scale boundaries), no saturation of the hierarchy 
< f (p) ^ J < f) and strong non universalities in 
the prefactors Aj m . 

We first present in Fig. 2 results concerning the isotropic 
sector, j = 0, m = 0, comparing the undecomposed struc- 
ture functions in the three direction with the projection 
S'p (|r|) and their logarithmic local slopes (inset). Only 
for the projected correlation it is possible to measure 
(5% of accuracy) the scaling exponents by a direct log- 
log fit versus the scale separation, \r\. The best fit gives 
£J=°(2) = 0.70±0.03. The undecomposed structure func- 
tions are overwhelmed by anisotropic effects at all scales 
which spoil completely the scaling behavior. 

In Fig. 3 we present an overview of all sectors j, m 
which have a signal-to-noise ratio high enough to ensure 
stable results [Q. Sectors with odd js are absent due 
to the parity symmetry of our observable. We measure 
anisotropic fluctuations up to j = 6. We notice from Fig. 
3 a clear foliation in terms of the j index: sectors with the 
same j but different ms behaves very similarly. In Table 
1 we present a more quantitative analysis by showing the 
results for the best power law fit for structure functions 
of orders p = 2,4. The first result we want to notice 
is the absence of any saturation for the exponents as a 
function of the j value. Unfortunatly the presence of an 
oscillation in all j ' = 2 sectors prevents us from measuring 
with accuracy the exponents in this sector, we therefore 
refrain from giving any number in this case. 

Let us also notice that the values for j = 4 and j = 6 
are different from what one would have expected if the 
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anisotropic effects would be given by simple smooth large 
scale fluctuations (see Table f). This fact leads to the 
conclusion that anisotropies are certainly the results of 
non-linear interactions in our flows, whether they corre- 
sponds to "homogeneous" fluctuations like in the Kraich- 
nan models or to some dimensional balancing between 
the non-linear terms and the forcing term is still an open 
question. The presence of a hierarchical monotonic in- 
creasing of exponents at fixing p and changing j leads 
to the possible breaking of the locality assumption in 
high enough j sectors For locality here we mean 

the fact that all integrals of pressure- velocity correlation 
functions are convergent both in the 1R and in the UV 
limits. 

Let us conclude by assessing also the important point 
connected to the existence of intermittency in higher j 
sectors. From Table 1 we see that already for the j = 4 
sector, and even more for j = 6, the fourth order scaling 
exponents are "almost" saturated, i.e. very close to the 
values of the second order exponents, ft is hard to say 
how much such a result is a quantitative sign of strong 
intermittency, due to the fact that we lack a clear un- 
ambiguous dimensional -non intermittent- prediction for 
anisotropic exponents (see below). A fast saturation of 
exponents within each sector as a function of the order of 
the moment must somehow be expected. We imagine the 
statistics in the anisotropic sectors being strongly domi- 
nated by "persistent" large scale structures, introducing 
cliffs structures (statistically speaking) characteristic also 
of saturation of exponents in anisotropic scalar advection 
p2Lp3l. Saturation, in the anisotropic sectors as a func- 
tion of the order of the observed moment, p, may also lead 
to the appearance of "persistency of anisotropies" even 
in presence of the observed strict hierarchical ordering, 
(&(p) < £ J (p) if J < j') as remarked in ||. In con- 
clusion we have presented the first numerical exploration 
of an anisotropic homogeneous turbulent flow. We have 
confirmed that by decomposing longitudinal structure 
functions in terms of the eigenvectors of the rotational 
operator we have a dramatic improvement of the scal- 
ing behavior in the isotropic sector. We have also used 
the SO (3) decomposition in order to asses two important 
questions opened in the field of anisotropic turbulence: 
(i) the presence of a hierarchical organization of turbu- 
lent fluctuations as a functions of the degree of anisotropy 
labeled by the j index (ii) the existence of intermittency 
(saturation as a function of the structure function order) 
in anisotropic sectors. The numerically measured values 
for the scaling exponents £ J (p) are not consistent with a 
simple "smooth" hypothesis for the nature of anisotropic 
fluctuations. More work is needed in order to understand 
the universality degree of our results as a function of the 
anisotropic properties of the large scale forcing. More 
work will be also devoted to measure fully tensorial quan- 
tities like Dij(r) = {(vi(x)—Vi(x + r))(vj(x) — Vj(x+r))) 
in order to be able to probe also odds sectors of the SO (3) 



group. 

We conclude by noticing that dimensional prediction 
for the £ 3 (p) with j > are far from being trivial. Indeed, 
different dimensionless quantities can be built by using 
some anisotropic mean observable (the mean shear for ex- 
ample, or the mean squared shear in our RKF) and the 
energy dissipation. Dimensional predictions than, would 
depend on the requirement that the anisotropic correc- 
tion is (or is not) an analytical, smooth deviation from 
the isotropic sector. 
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(j, m ) 


(0,0) 


(4,0) 


(4,2) 


(6,0) 


(6,2) 


61*2 

€4|Cf 


0.70 ± 0.03 | 2 
1.28 ± 0.05 | 4 


1.67 ± 0.07 | 2 
2.15 ± 0.1 | 4 


1.7 ±0.1 | 2 
2.2 ±0.1 | 4 


3.4 ± 0.2 | 4 
3.2 ± 0.2 4 


3.3 ± 0.2 | 4 
3.2 ± 0.2 4 



TABLE I. Best fit of the scaling exponents in all stable sectors. For comparison we also give, £p , the exponents for the 
case of a smooth (many times differentiable) anisotropic field. Some sectors are absent due either to the small signal-to-noise 
ratio or to the presence of sign changes in S^ m (|r|) which prevent the very definition of a slope. Errors are estimates from the 
fluctuation of the logarithmic local slopes at resolution 256 3 . 
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FIG. 1. Typical energy (above) and energy dissipation 

(below) time evolution in arbitrary units of the Ran- 

dom-Kolmogorov flow at resolution L x = L y — L z = 256. 

Inset: root mean squared velocity (v x ) as a function of 

the spatial location in the three directions: (v x (x/L x )) (x), 

( V x(y/L y )) (*), (v x (z/L z )) (+). For comparison is also shown 

the same quantity (□) from experimental state-of-the-art 

anisotropic homogeneous shear flow at changing the position 

along the shear direction y [2]. All curves are normalized to 

be 1 at x/L x = y/L y = z/L z — 0.5. 
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FIG. 2. Isotropic sector. Log-log plot of ^'"(M) versus 
\r\ (+), and the three undecomposed longitudinal structure 
functions in the three directions x,y,z (d,*,x) respectively, 
at resolution 256 3 . The straight line has the best fit slope 
£-' =0 (2) _ o 7Q Inset: logarithmic local slopes of all curves 
(same symbols) plus the straight line corresponding to the 
intermittent isotropic scaling. Notice the dramatic improve- 
ment in the scaling behavior of the projected correlation. Sim- 
ilar results hold for higher orders p > 2 (not shown). 
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FIG. 3. Log-log plot for the absolute value of the projected 

second order structure functions, |>S2' m (M)|, vers us the scale 
r, on all measurable sectors (up to j = 6). Sectors: (0, 0), (+); 
(2,2), (X); (4,0), (□); (4,2), (*); (6,0), (o); (6,2), (■). The 
statistical and numerical noise induced by the SO(3) projec- 
tion can be estimate as the threshold where the j = 6 sector 
starts to deviate from the monotonic decreasing behavior, i.e. 
O(10" 3 ). 
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